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Abstract 

It is shown that the operator methods of super symmetric quantum mechanics and the 
concept of shape invariance can profitably be used to derive properties of spherical harmon- 
ics in a simple way. The same operator techniques can also be applied to several problems 
with non-central vector and scalar potentials. As examples, we analyze the bound state 
spectra of an electron in a Coulomb plus an Aharonov-Bohm field and/or in the magnetic 
field of a Dirac monopole. 
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Spherical harmonics are introduced in physics courses while treating the Laplacian in 
spherical polar coordinates. In the context of quantum mechanics, the Schrodinger equa- 
tion is separable into radial and angular parts if the potential is spherically symmetric. 
The angular piece of the Laplacian operator generates spherical harmonics which obey 
interesting and useful recursive formulae. With the advent of supersymmetric quantum 
mechanics (SUSYQM)[l-3] and the idea of shape invar iance[U, study of potential prob- 
lems in nonrelativistic quantum theory has received renewed interest. SUSYQM allows one 
to determine eigenstates of known analytically solvable potentials using algebraic operator 
formalism Q without ever having to solve the Schrodinger differential equation by stan- 
dard series method. However, the operator method has so far been applied only to one 
dimensional and spherically symmetric three dimensional problems. 

The study of exact solutions of the Schrodinger equation with a vector potential and a 
non-central scalar potential is of considerable interest. In recent years, numerous studies [6- 
10] have been made in analyzing the bound states of an electron in a Coulomb field 
with simultaneous presence of Aharonov-Bohm (AB) [Tl|l field and/or a magnetic Dirac 
monopole|l2|]. In most of these studies, the eigenvalues and eigenfunctions are obtained via 
separation of variables in spherical or other orthogonal curvilinear coordinate systems. The 
purpose of this paper is to illustrate that the idea of supersymmetry and shape invariance 
can be used to obtain exact solutions of such non-central but separable potentials in an 
algebraic fashion. In this method, it emerges that the angular part (as well as the radial 
part) of the Laplacian of the Schrodinger equation can indeed be dealt with using the idea 
of shape invariance. This is a novel method of generating interesting recurrence properties 
of spherical harmonics. In standard text books on quantum mechanics^], properties of 
spherical harmonics are discussed from a different point of view. In this regard our ap- 
proach is new and instructive in the sense that the radial and the angular pieces of the 
Schrodinger equation can both be treated within the same framework. The present work 
also gives insight into the solvability of certain three-dimensional problems. Basically, this 
is a consequence of both separability of variables and the shape invariance of the resulting 
one-dimensional problems. 
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The Schrodinger equation for a particle of charge e in the presence of a scalar potential 
V(r, 9, 4>) and a vector potential A(r, 9, (f>) is (in units of H = 2m = 1) 



[{-iV - eA(r, 9, <j))} 2 + V(r, 9, <f>)]i(> = Eifj 



F(0) 

We will consider vector potentials of the form A = -4^4 , and scalar potentials of the 
form V = V±(r) + ^ ^ ■ The Schrodinger equation reads 

1 d 2 ^ 
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r 2 sin 2 i 



+ V 1 {r) + 



V 2 {9)\ 2ieF(9)dip 



r 2 sin 2 i 



(2) 



Eq. @ permits a solution via separation of variables, if one writes the wave function in the 
form 

i/>(r,6,4>) = R(r)P(9)e imrp . (3) 
The equations satisfied by the quantities R(r) and P{9) are: 



d 2 R | 2dR | / ^ /(/ + !) ' ;> () 

dr 2 r dr \ ' ' r 2 



d 2 P dP 



sin 



P = 0. 



(4) 



(5) 



Solutions of these equations for various choices of scalar and vector potentials will now be 
discussed. 



Spherical Harmonics: Let us begin with the simplest case of a free particle [Vi(r) = 
V2{9) = F(9) = 0.] This will allow us to obtain the standard properties of the spherical 
harmonics. In this case, differential equation (]5|) reduces to 

,2 



d 2 P dP 



1(1 + 1) 



m 



sin 



P = 



(6) 



which is the equation satisfied by associated Legendre polynomials. To solve it by the 
SUSYQM method, we need to re-cast it into a Schrodinger-like equation. Changing the 
variable 9 — » z through a mapping function 9 = f(z), one obtains 



£P 

dz 2 



+ 



/ 



7+ /'cot/ 



dP 

dz 



+ f 
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1(1 + 1) 



m 



sin 2 / 



P = . 



(7) 
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Since we want the first derivative to vanish, we choose 

"/'cot/ , (8) 

which gives 

9 = f = 2tan~ 1 (e 2 ) . (9) 

This transformation amounts to the replacement sin9 = sechz and cos 9 = — tanhz. The 
range of the variable z is — oo < z < oo. Eq. now reads 

d 2 P 

--To -1(1 + l)sech 2 z P = -m 2 P . (10) 

dz A 

This is a well-known, shape invariant, exactly solvable potential. It can be readily solved 
by standard SUSYQM techniques |, g. The energy eigenvalues for the potential V(z) = 
—ao(a,Q + l)sech 2 z (ao > 0) are: 

E n = -(a -n) 2 ; (n = 0, 1, 2, • ■ ■ JV) , (11) 

where N is the number of bound states this potential holds, and is equal to the largest 
integer contained in ao- The eigenfunctions Tp n (z, oq) are obtained by using supersymmetry 
operators ||]: 

tp n (z;a ) oc A\z;a )A\z;a 1 ) ■ ■ ■ A 1 (z; a n -i)ifj (z; a n ) , (12) 

where A\z;a) = ( — ^ + atanhz^ and a n = ao — n. The ground state wave function is 
ipo(z; a n ) = sech a "2. 

For our problem, ao = I, E n = —m 2 and consequently one has 

n = l-m; Pj j7n (tanh.z) ~ ^i- m (z; I), (13) 

where P^ m (tanh z), the solutions of eq. (|6|), are the associated Legendre polynomials of 
degree I. Now that these Pi m (z) functions can be viewed as solutions of a Schrodinger 
equation, we can apply all the machinery one uses for a quantum mechanical problem. 
For example, we know that the parity of the n-th eigenfunction of a symmetric potential 
is given by (—1)™, we readily deduce the parity of P^ m to be (— \) l ^ m = (— l)' +m p5[|. 
Also, the application of supersymmetry algebra results in identities that are either not very 
well known or not easily available. With repeated application of the A* operators, we can 



determine P^ m for a fixed value (I — m). As an illustration, we explicitly work out all the 
polynomials for I — m = 2. The lowest polynomial corresponds to I = 2, m = 0. For a 
general I, using eq. (|i~2]) , one gets 

■ A t (5f;0A t («;Z-l)^(2;i-2) 



P/ ) /_ 2 (tanhz) 



d 



dz 



d 



- — h Z tanh z I —3 — h (Z — 1) tanh z sech 



dz 



J-2. 



-1 + (2/ - l)tanh 2 z 



sech 



Z-2. 



(14) 



A similar procedure is readily applicable for other values of n. When one converts back to 
the original variable 9 (sechz = sin# ;tanhz = — cos#), the results are: 



-Pz.z (cos ( 
Pz,z-i(cos( 

PZ,Z-2(C0S( 
Pz,Z-3(cOSl 
-Pz,Z-4(cOS( 

P/ ; z_ 5 (cos! 



sin 



~ sin' ^ 9 cos 1 



(15) 



1-2, 



sin 



sin' 3 9 cos 1 



-1 + (21 - l)cos 2 6> 
-3 + (2/ - 1) cos 2 9 
3 - 6(2/ - 3) cos 2 9 + (21- 3)(2Z - 1) cos 4 9 
15 - 10(2/ - 3) cos 2 9 + (21- 3) (2/ - 1) cos 4 



sin 



Z-4, 



sin 



Z-5, 



COS I 



These results are not readily available in standard books, although of course they can be 
obtained from the generating function for associated Legendre polynomials. Likewise, there 
are several recurrence relations which are easily obtainable via SUSYQM methods. In 
particular, by applying A or A] once, we generate recurrence relations of varying degrees 
(differ in /): 

d 



ipi- m (z,l) = Ai(z,l)ipi- m -i(z,l - 1) 
Tpi- m (z, I) = A(z, I - l)^_ m+1 (z, / + 1) 



Pl,m(x) 



'l-x')-^ + lx)P^ 1>m (16) 



P Lm (.r) I (l-x 2 )^ + (/-l)x| AV,,„. 



Non-Zero Vector Potential: Here we examine the general case with non-zero F(9) in 
eq. (||). In particular, we consider 



F(9) = ^- + g(l-cos9) , V 1 (r) = -— , V 2 (9) = 
2tt r 



(17) 



This choice corresponds to the physically interesting problem of the motion of an electron 
in a Coulomb field in the presence of an Aharonov-Bohm potential Aab = 2irrsm0 ^<t» an< ^ 
a Dirac monopole potential Ajj = S^e^jlC]. Both potentials can be chosen to be of 



arbitrary strength depending on the values of the coupling constants T and g. Again, we 
transform eq. (|5|) via the change of variables from 9 to z given in eq. @. The result is 

d 2 P 



dz 2 



+ [(X 2 - rh 2 ) - V{z)\P = 



with 



and 



V(z) = (A 2 + q 2 ) tanh 2 z — 2qm tanh z 



-ge , m 



2^ 



+ ge-m , X 2 = 1(1 + 1) 



(18) 
(19) 
(20) 



The potential V(z) has the form of the Rosen- Morse II potential, which is well-known to 
be shape invariant. More specifically, the potential 

Vrm(z) = a o( a o + 1) tanh 2 z + 2b tanh z (6 < o,q) (21) 



has energy eigenvalues and eigenfunctionsl 



IpT) 



a (a + 1) - (a - nf 



h 2 



(1 - tanhz) : 



ao-n+ 



(a - n) 2 ' 



[1 + tanh z) 



(22) 



I ao-nH 



r,ao 



a -n J 



For our case, E n = X 2 — rh 2 and the constants ao, and bo are given by 

a 



-\ + ^\ + (X 2 + q 2 ) 1 



-qm 



Using these values and stipulating that eigenfunctions given in eq. ( |22D be normalizable, 
we get 



n 



\ + (X 2 + q 2 ) 



Corresponding to these two cases, one gets 

i 



m 



I 



+ \ (n+ q + 



1^2 



+ J(n+\fh\ + \) 2 



(\q\ > \m\) 
(\q\ < |m|) 

(|g| > |m|) 
(\q\ < \rh\) 



(23) 



(24) 



The energy eigenvalues obtained from eq. 

En 



for the Coulomb potential are 



(25) 



A[N + l + l] 2 ' 

Therefore, our final eigenvalues for a bound electron in a Coulomb potential as well a 
combination of Aharonov-Bohm and Dirac monopole vector potentials are 

Z 2 e 4 



En 



-Z 2 e A 



JV +i+v / ( n +H+|) 2 -9 2 



(M > H) 

(|g| < |m|) 



(26) 



which agree with eqs. (32) and (33) respectively of ref. (K 

The above calculation using operator techniques in SUSYQM can also be carried out 
when the Coulomb potential is replaced by a harmonic oscillator. Potentials of this type 
have been studied by many authors using other approaches ||. 

Non-Central Scalar Potential: As a final example, we consider the case of zero vector 
potential and a scalar piece consisting of the Coulomb potential and a non-central part V2(0). 
The angular piece satisfies a modified version of eq. ([l0|): 
d 2 P 



dz 2 



+ 



sech^ {1(1 + 1) - V 2 (z)} - mf 



P = 



(27) 



This equation is the Schrodinger equation for one of the known shape invariant potentials 
provided the function V2(z) is chosen appropriately. The following three simple choices 
can be made: V2(z) = ftsinhz, &sinh2z, 6cosh 2 z. These choices give rise to the Scarf II 
and Rosen-Morse II potentials)^] in eq. (27). In terms of 9, these choices correspond to 
non-central potentials with angular dependences cot 6, cot #cosec# and cosec 2 #, respectively. 

Here, we will treat the case V^(^) = bsinhz in detail. The full potential in eq. ( p7|) is 
— 1(1 + l)sech 2 z + teechztanhz and the role of energy is played by —m 2 . RecallQ that the 
potential V (x) = (B 2 — A 2 — yl)sech 2 x + B(2A + l)sechx tanh x with A > has eigenvalues 
E n = 2 An — n 2 (n < A) . In terms of our parameters, this implies 



£(2,4 + 1) = 6 , 
Elimination of A, B leads to 



b 2 -A 2 



A 



-1(1 + 1) 



2An - n 



-m 



(28) 
(29) 



with 

(n 2 — m?) (n 2 — m 2 ) 2 n 2 b 2 

2n An 2 (n 2 — m 2 + n) 2 

Substitution for I into eq. ( p5[ ) gives the eigenvalues for this problem. Although we have ex- 
plicit energy eigenvalues and eigenfunctions, it is not clear whether the non-central potential 
we have considered has physical significance. 

In conclusion, we mention that in this paper we have attempted to explore the effective- 
ness of the SUSYQM operator methods to obtain analytic solutions of Schrodinger systems 
in more than one dimensions. It is clear that the factorization technique which has so 
far been applied to only one dimensional or spherically symmetric three dimensional prob- 
lems, can be equally useful for non-central separable problems for which one dimensional 
equations can be recast into Schrodinger equations with shape invariant potentials. Many 
interesting properties of spherical harmonics emerge naturally as a simple realization of this 
operator technique. One should also note that although, in this paper, we have focused on 
spherical polar coordinates, any orthogonal curvilinear coordinate system which is separable 
into Schrodinger-type equations with shape invariant potentials will allow similar algebraic 
analysis, and will have analytically solvable eigenvalues. 
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